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SUMMARY 

We consider classical estimators for a class of physically realizable linear quantum systems. Optimal 
estimation using a complex Kalman filter for this problem has been previously explored. Here, we study 
robust Hoo estimation for uncertain linear quantum systems. The estimation problem is solved by converting 
it to a suitably scaled Hoo control problem. The solution is obtained in the form of two algebraic Riccati 
equations. Relevant examples involving dynamic squeezers are presented to illustrate the efficacy of our 
method. Copyright © 2016 John Wiley & Sons, Ltd. 
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1. INTRODUCTION 

In recent years, there has been significant interest in studying estimation and control problems for 
quantum systems [1-13]. Linear quantum systems are an important class of quantum systems and 
have been of particular interest in this context [1,4—7, 12,14-20]. Such linear quantum systems 
have been useful in describing quantum optical devices such as linear quantum amplifiers [15], 
finite bandwidth squeezers [15] and optical cavities [14,21]. Coherent feedback control has also 
been studied a lot recently for linear quantum systems [5, 6,8,9,18,20, 22-24]. The authors have 
previously explored a related coherent-classical estimation problem [25,26], where the estimator 
consists of a classical part, that produces the desired final estimate and a quantum part, which may 
also involve coherent feedback. In this work, the authors have studied optimal classical estimation 
for linear quantum systems using a complex Kalman filter. 

Ref. [27] considered a quantum observer, that is a purely quantum system, which produces a 
quantum estimate of a variable for the quantum plant. By contrast, we here consider classical 
estimation for linear quantum systems, where the estimator is a classical system that yields a 
classical estimate of a variable for the quantum plant. A robust quantum observer for uncertain 
quantum systems was constructed in Ref. [4]. On the other hand, here we build a robust complex 
classical estimator for uncertain linear quantum systems. A robust classical Hoo estimator for 
uncertain linear systems was presented in Ref. [28]. Here, we develop a more general and complex 
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Hoc filter for robust classical estimation of uncertain linear quantum systems. The solution to the 
Hoc estimation problem is obtained by means of two algebraic Riccati equations, upon converting 
the estimation problem to a scaled Hoc control problem. 

The paper is structured as follows. We introduce the class of linear quantum systems considered 
in this paper in Section 2. The complex Kalman filter is discussed in Section 3 for optimal classical 
estimation of linear quantum systems. Section 4 then considers the robust Hoc estimation problem 
for uncertain linear quantum systems and presents our main result. Illustrative examples of our 
robust estimator for two different scenarios are provided in Sections 5 and 6. Finally, Section 7 
concludes the paper with relevant remarks and possible future work. 


2. LINEAR QUANTUM SYSTEMS 


The class of linear quantum systems we consider here are described by the quantum stochastic 
differential equations (QSDEs) [25]: 


where 


da(t) 

da{t)* 

dy(t ) 
dy(t)# 


= A 

i i 

^ a 

C-+* ^-s 

dt + B 

i i 

,-.=8= 

If 

= C 

a(t) 

dt + D 

dA(t) 

dA{t)* 


A=A(A ll A 2 ), B = A{B u B 2 ), 
C = A(C 1 ,C 2 ), D = A{D u D 2 ). 


( 1 ) 


( 2 ) 


Here, a(t) = [a\ (t)... a n (t)] T is a vector of annihilation operators. The adjoint of the operator 
<H is called a creation operator, denoted by a*. The vector A(t) = [A\ it)... A, n (t)\ r represents 
a collection of external independent quantum field operators and the vector y represents the 
corresponding vector of output field operators. Also, the notation A(Ai, A 2 ) denotes the matrix 


. Here, A u A 2 £ C nxn , B 1 , B 2 £ C nxm , C u C 2 £ C mxn , and D i, D 2 £ 


A\ A 2 
A* A* 

Moreover, # denotes the adjoint of a vector of operators or the complex conjugate of a complex 
matrix. Furthermore, t denotes the adjoint transpose of a vector of operators or the complex 
conjugate transpose of a complex matrix. 


Definition 2.1. [25] A complex linear quantum system of the form (1), (2) is said to be physically 
realizable if there exists a complex commutation matrix 0 = 0^, a complex Hamiltonian matrix 
M = M ', and a coupling matrix N such that 


0 = tjt\ 


(3) 


where J = 


/ 0 
0 -I 


,T = A(Tj, T 2 ) is non-singular, M and N are of the form 


M = A(M 1 ,M 2 ), N = A{N u N 2 ), 


(4) 


A = —lQM — 

b = -erfj, 
C = N, 

D = I. 


(5) 
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Here, t = yf-A is the imaginary unit, and the commutation matrix 0 satisfies the following 
commutation relation: 


a 

_ a* \ ’ 

a 

a* _ 

a 

~ a* 

a 

a# 


= ©. 



( 6 ) 


One can verify that 0 is a 2n x 2n matrix, the elements of which are as follows, given 
i,j = l...n: 


0 'ij — [Vx,, (i.j ], 

0(n+i)(n+ i j) = 

®i(n+j) = 

0 (n+l)i = [ala*], 
= 1, 

®(n+X)(n+i) = 


The annihilation and creation operators can be used to construct the number operator N = a^a, 
the eigenstates of which form the orthonormal number (or Fock) states [21]: 

N\q) = cJa\q) = q\q), q = 0,1,2,... 


In particular, the state |0) is called the vacuum state [21]: 


a|0) = 0. 

The annihilation and creation operators have the properties of lowering and raising the number 
of a state respectively [21]: 


a\q) = Vq\q- 1), 
a*\q) = y/q+l\q + l). 


Theorem 2.1 

[25] The linear quantum system (1), (2) is physically realizable if and only if there exists a complex 
matrix 0 = 0i such that 0 is of the form in (3), and 

AQ + 0x4 f + BJB^ = 0, 

B = -ecu, (7) 

D = 1. 

If the system (1) is physically realizable, then the matrices M and N define a complex open 
harmonic oscillator with a Hamiltonian operator 

and a coupling operator 

L = [ iVi 


a T ] M 


N 2 


a 

a# 
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Figure 1. Schematic diagram of classical estimation for a quantum plant. 


3. KALMAN FILTER 

The schematic diagram of a classical estimation scheme is provided in Fig. 1. We consider a 
quantum plant, which is a quantum system of the form (1), (2), defined as follows: 


da(t) 

da{t)* 

= A 

a{t ) 

_ a{t)* 

dt. + B 

■ i 

If 

dy(t) 

_ dy(t)# 

= C 

a{t) 

. 

dt + D 

dA(t) 

dAit)* 


r = L 

a(t) 
a{t)* _ 




Here, z denotes a scalar operator on the underlying Hilbert space and represents the quantity to 
be estimated. Also, y(t) represents the vector of output fields of the plant, and A(t) represents a 
vector of quantum disturbances acting on the plant. 

A quadrature of each component of the vector y(t) is measured using homodyne detection to 
produce a corresponding classical signal yp. 


dyi 


e -^i 

~VT 


piOi 

dy 1 + -^dyi 


dy m 


V2 dym + 



(9) 


Here, the angles &i,9 m determine the quadrature measured by each homodyne detector. The 
vector of classical signals y = [y-\ ... y m ] T is then used as the input to a classical estimator defined 
as follows: 


dx e = A e x e dt + K e dy, 

e e (10) 

z = L e x e . 

For the sake of comparison, we will first consider the optimal estimation problem for quantum 
linear systems; see also [4,29]. The optimal classical estimator is given by the standard (complex) 
Kalman filter defined for the system (8), (9). This optimal classical estimator is obtained from the 
solution to the algebraic Riccati equation: 

AP + PA* + BB* - {B + PC^S^SiB + PC t ) t = 0, (11) 
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where 


S = [ Si S 2 ] , 


Si 




s 2 


V2 

0 



0 

0 


j 

0 

0 



( 12 ) 


Here we have assumed that the quantum disturbance A is purely canonical, i.e. dAdJ\) = Idt 
and hence D = I. 

Then, the corresponding optimal classical estimator (10) is defined by the equations: 


A e = A - K e SC, 

K e = {B + PC^)S\ (13) 

L e = L. 


4. ROBUST Hao FILTER 


Corresponding to the system described by (8), (9), we define our uncertain system modelled as 
follows: 

(Si) : x(t) = [A + AA(t)]x(t) + [B + AB(t)]w(t), 

z(t) = Lx(t), (14) 

y\t) = S[C + AC(t)]x(t) + SDw(t ), 


where x(t) := 


a(f) 

a(t)# 


is the state, w(t) is the disturbance input, z(t) is a linear combination 


of the state variables to be estimated, y'{t) is the measured output, L £ O ' x2n , SC £ C'" x2n , 
SD £ C mx2m , and A,4(-), A B(-) and A C(-) denote the time-varying parameter uncertainties. 
These uncertainties are in the following structure 


' A A(t) ' 


' H x ' 

AC(t) 


. H z . 




AB(t) = H 2 F 2 (t)G, 


05) 


where Hi, H 2 , H 3 , E and G are known complex constant matrices with appropriate dimensions, 
and the unknown matrix functions Fi(-) and F 2 (-) satisfy the following: 

F}(t)Fi{t)< I, Vi, 

+ (lu) 

F}(t)F 2 (t)<I, Vi. 


Note that for the system (14) to be physically realizable, the following constraints are required to 
be satisfied by the uncertainties: 


AA0 + 0AA f + BJAB ] + A BJB X + ABJAB^ = 0, 

A B = -0A C ] J. 


(17) 
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The robust // lOC estimation problem for the uncertain system ( 14) can be converted into a scaled 
Hoc control problem, as described in Ref. [28], by introducing the following parameterized linear 
time-invariant system corresponding to system (14): 

(S 2 ) : x(t) = Ax(t ) + [ B ^Hi A-H 2 ] w(t), 



ei E 


1 

O 

O 

O 

_1 


0 

5(f) = 

0 

L 

x(t) + 

1 

Ct\ 

O O 

O O 

1 _ 

w(t) + 

0 

-I 


y'(t ) = SCx(t) + [ SD £ SH 3 0 ] w(t). 

Here, u(t) is the control input, 5(f) is the controlled output, ei, e 2 > 0 are suitably chosen scaling 
parameters and 7 > 0 is the desired level of disturbance attenuation for the robust H aa estimation 
problem. We also have the augmented disturbance 


w(t) := 


w(t) 


where 

77(f) := F 1 (t)Ex{t ), 

and 

m ■= F 2 (t)Gw(t). 

The following assumptions are made for the system (18): 


Assumptions 

Al. The system matrix A is stable. 

A2. e^G+GcJ. 

A3. [ SD SH 3 ] is of full row rank. 


A4. rank 


A — lujI B 
C D 


= 2 n + 2m, 


Vw > 0. 


Remark 1. The assumption Al is required so that the Hoc norm for the combined plant-estimator 
system is finite. The remaining assumptions are technical assumptions arising in Hoc control theory 
which are required in order to obtain a solution using Riccati equations. 

A complete solution to the robust Hoc estimation problem is then provided below. 

Theorem 4.1 

Consider the robust Hoc estimation problem for the uncertain system (14) converted to a scaled I 
control problem for the system (18), satisfying the assumptions Al to A4. Given a prescribed level 
of disturbance attenuation 7 > 0, the robust Hoc estimation problem for the uncertain system (14) 
is solvable if for some n, e 2 > 0 , the following conditions are satisfied: 

(a) There exists a stabilising solution X = > 0 to the algebraic Riccati equation: 


Al X + XA + A( 7 " 2 BiB 1 )A + c\(I - D 12 E 1 1 d\ 2 )C 1 = 0 . (19) 


(b) There exists a stabilising solution Y = > 0 to the algebraic Riccati equation: 

AY + YA f + YclCiY + 7" 2 BiB t 1 

- ( 7 - 1 b 1 ^ 1 + 1 ycI)s% 1 s( 1 ~ 1 b 1 d1 1 + 7 y^)t = 0 . 


( 20 ) 


(c) I - 7 ~ 2 XY > 0. 
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Here, we have 


A = A, 

Hi = [ B(I - elG^G)- 1 / 2 


C 1 = 

D\2 = 


e\E 

0 

L 

0 

0 

-I 


c 2 = c, 


D 21 =[ Dil-elG^G)- 1 ' 2 
S=S, 


l-Ih 

ei L 


2-h 2 

C2 Z 


;H 3 0 ] 


and 

Ei = d\ 2 Di2 = I , 

2 

E 2 = ~SD 2 iD\J? = SD(I - e^G t G)" 1 L> t S' t + ^-S , H 3 H 3 t 5' t . 

6 1 

When conditions (a)-(c) are satisfied, a suitable estimator is given by: 


x{t) = A K x(t ) + B K y'(t ), 
z(t) = C K x(t), 


where 

Art =A~ BkSC 2 + 7 " 2 (^i - B k SD 2 i)b\x, 
b k = 7 2 (/ - + 7 - 2 H 1 ^ 1 H t )H“ 1 , 

Gk = 

Proof 

The system (18) is of the following form: 

(E 2 ) : x(t) = Ax(t) + Biw(t) + B 2 u(t), 

z(t) = Cix(t) + Duw{t) + Di 2 u(t), 
y'(t) = SC2x(t) + SD 2 iw(t ) + SD 22 u(t), 


where 


Bi=[B 1 -Hi %H 2 ], 

b 2 = 0, 


Cl = 

Hn = 

D\2 = 


Cl E 
0 , 

L 

0 0 0 

e 2 G 0 0 

0 0 0 

0 

0 , 

I 


c 2 = c, 

D 2 1=[D ^H 3 0 ] , 

D 22 = 0. 


( 21 ) 


( 22 ) 


(23) 


(24) 


(25) 


(26) 
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We will use the results from Ref. [30] to solve the above control problem. However, Ref. [30] 
requires the matrix Du to be zero. The system (18) with non-zero Du can be converted to an 
equivalent system having Du = 0 using the loop-shifting technique, as outlined in Ref. [31]. 

We define 


Dun D m 2 

D\l2l Du22 + Dqq ’ 


(27) 


where we let D x = -Du 22 - £>1121 (7 - D\ ni Diin) 1 nt 


we take Dun = 


0 

e 2 G 


1111-^1111 
, Du21 = [ 0 0], Z?1112 = 


Z?] m i7iii2. Here, for our system (18), 


and D 1122 = 0. Note that 


| D 11 11 < 1 follows from Assumption A2. Then, we get Doe = 0. Thus, we have Du = D\ 
Hence, the //^ control problem in (25) takes the following form: 


(S 3 ) : x(t) = Ax{t ) + Biui(t) + B 2 u{t), 

z(t) = C\x{t) + Duw(t) + D 12 u(t), (28) 

y'[t) = SC 2 x(t) + SD 21 w(t) + SD 22 u(t), 


where 

A = A B 2 D 00 C 2 = A, 

B\ = B\ + B 2 D 00 D 2 \ = Bi, 

B 2 = B 2 = 0, 

Ci = Ci + Di 2 D 00 C 2 = Ci, 

Di 2 = D 12 , (29) 

C 2 = C 2 , 

D 2 i = D 2 i, 

D 22 = D 22 = 0, 

S = S. 

The Hoo control problem equivalent to the above system (28) is [31]: 

(E4) : x[t) = Ax(t) + Biw[t) + B 2 u{t), 

z(t) = Cix(t) + Duw{t) + Di 2 u(t), (30) 

y'(t) = SC 2 x(t) + SD 2 iw(t) + SD 22 u{t), 


Here, 


A = A + BiR^D^Ci, 

_ D TD~ 1/2 

&l — iD\ri x , 

B 2 = B 2 + BiR 1 l L)\ 1 Di 2 , 
Ci = r~ 1/2 Ci, 

Du = 0 , 

D u = Rf 1,2 Di 2 , 

c 2 = c 2 + d 2 iR 1 1 b\ 1 Ci, 

b 2 i = b 2 iR- 1/2 , 

D 22 = £>217?! 1 l)} 1 Di2, 

S = S, 


(31) 
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where 


R\ — I — -D11-D11 — 

Ri = I- -D11-D11 = 


I — e\G^G 0 
0 I 

0 0 


I 0 

0 I- e|GG t 
0 0 


0 

0 

/ 

0 

0 

I 


(32) 


One can verify that (31) yields (21), B 2 = 0 and D 2 2 = 0. Note that from assumptions A1 to A4 
and (22), it follows that we have: 






Ei > 0 . 


E 2 > 0. 


rank 

rank 


A — lojI 

b 2 

Ci 

D12 

A — uyl 

Bi 

c 2 

D21 


= 2 n + p for all u > 0. 


= 2 n + m for all ui > 0. 


Hence, the control problem for the system (30) can be solved using the results in Ref. [30]. 
However, Ref. [30] assumed 7 = 1 . The results from that paper may be generalised for different 
values of 7 > 0, simply by scaling the coefficients of the disturbance ui(t) in (30), viz. B 1 and D 2 1 
(note .Du = 0), by 7 -1 . Note that this also has an effect on E 2 , which is scaled by 7 -2 . This yields 
conditions (a)-(c) of the theorem, which are required to be satisfied by the system (30), such that a 
suitable estimator is given by (23). □ 


The transfer function of the robust estimator (23) can be obtained to be: 


Gk(s):=4A=C k (sI-A k )- 1 B k . (33) 

y (s) 

The estimation error is given as: 

e(t) := z(t) — z(t) = Cxx{t) — Lx(t). (34) 

Then, the disturbance-to-error transfer function may be obtained as: 

G we (s) := = [ -L C K } (al- 

ix {s) L J \ 

We are interested in the disturbance A to error e transfer function, which is simply the first 
component of the matrix transfer function G we (s). The other component is the disturbance to 
error e transfer function, which we shall ignore. 

Remark 2. Note that the fact that the plant is a quantum system that will be physically realizable 
restricts the class of plants under consideration, when compared to the case when the plant is a 
classical system as in Ref. [28], owing to the conditions in (7) (and also (17)) required to be satisfied 
by the system matrices of the quantum plant (uncertain plant). 


A + AA 0 

V 

B + AB 

B k S(C + AC) a k 

J 

B k SD 


(35) 


5. NUMERICAL EXAMPLE 1 

An example of a linear quantum system from quantum optics is a linearized dynamic squeezer. 
This corresponds to an optical cavity with a non-linear active medium inside. Let us consider the 
quantum plant to be a linearized dynamic squeezer, described by the QSDEs [26]: 
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da(t ) 
da(t)* 



a(t) 

a(t)* 


dt — \[k 


dA(t) 
dA{t)* ’ 


dy(t) 

dy(ty 



dt -f- 


dA(t) 

dA(t)* 


z(t) = [ 0.1 -0.1 


a(t) 

a(t)* 


(36) 


where /3 > 0 is the overall cavity loss, k > 0 determines the loss arising from the cavity mirrors, 
X € C quantifies the size of the non-linearity of the active medium, and a is a single annihilation 
operator of the cavity mode. 

Here, we choose (3 = 4, k = 4, and x = — 0 . 5 . Then, the above quantum system is physically 
realizable, since we have /3 = k. Moreover, we fix the homodyne detection angle at 90 °. We, thus, 
have the following: 


A = 

D = 


-2 

0.5 

1 0 
0 1 




0 

-2 


-0.1 


C = 


2 0 
0 2 




(37) 


We introduce uncertainty in the parameter a := y/R as follows: a a + /.«5a:, where <5| < 1 is 
an uncertain parameter and fj, € [0,1) determines the level of uncertainty. Then, we will have the 
following: 


AA = 

pda 2 ^ 

S 2 od 

2 

0 

A B = 

1 

On 

P 

O 

o 


pSa 2 

o 

~E 

o 

» 


A G = 

pSa 0 

0 pSa 




Then, we define the relevant matrices in (15) as follows: 


Fi(t) = 


F 2 (t) = 


E = 

G = 
H\ = 
H 2 = 
h 3 = 


S 0 0 0 ' 

0 (5 0 0 

0 0 6 2 0 

0 0 0 ( 5 2 

5 0 ' 

0 <5 ’ 


1 0 
0 1 ’ 


2 [MX 2 

0 

p 2 a 2 

0 

0 

2 pa 2 

0 p 2 a 2 

—pa 

0 



0 

—pa 

5 


—2 pa 

0 

O 

O 



0 —2 iia 0 0 


(38) 


(39) 
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One can then verify that we have A A = II\ F\ (t)E, A B = H 2 F 2 (t)G and AC = H :i F\ (t)E, as 
required in (15). We set the uncertainty level to /r = 0.1. Moreover, in our simulations, we choose a 
fixed value of 5 = 1. 

We now solve the associated H 00 estimation problem using the Riccati equation approach 
described in the previous section. We choose the desired disturbance attenuation level to be 
7 = 0.65. Then, the scaling parameters are suitably chosen to be ei = 0.2 and e 2 = 0.6. A robust 
Hoo estimator is obtained as in (23) with the following parameters: 


A k = 


0.1905 -1.4676 

-1.4676 0.1905 ’ 


B k 


-1.4717t 
1.47174 ’ 


C K = [ 0.1 -0.1 ] . 


(40) 


The transfer function (33) of the estimator is obtained to be: 

^ —0.2943ts- 0.37596 

K ^ ~ s 2 -0.381s-2.118' 


(41) 


Fig. 2 shows the disturbance to error transfer function (35) of the above robust biter. It also shows 
those of the optimal //^ biter and the Kalman biter for the uncertain system (14) with 5 = 1 in our 
example for comparison, which corresponds to the uncertain parameter taking on its maximum 
value. The optimal Kalman and Hoo biters are built for the nominal system of (14). 

Clearly, the robust H a0 biter provides for preferable disturbance attenuation and performance 
than the standard optimal biter or the Kalman biter for this value of the uncertain parameter S. 


6 . NUMERICAL EXAMPLE 2 


In this section, we consider another numerical example involving the dynamic squeezer plant (36), 
(37). Here, we do not consider uncertainty in (3 or k. Instead, we introduce uncertainty in the 
squeezing parameter x as follows: \ ~> X + /A’x, where again 4| < 1 is an uncertain parameter 


Bode Diagram 
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and 0 < n < 1 is the extent of uncertainty. Then, we will have the following: 


AA = 


A B = 


A C = 


0 ~n5x 

-fi5x 0 J ’ 

0 0 ' 

0 0 J ’ 

0 0 ' 

0 0 ■ 


(42) 


Note that the results from Ref. [28] are enough to construct the robust estimator for this scenario. 
When A B = 0 as in this case, has no impact on the estimator. 

Here, the relevant matrices may be defined as follows: 


Fi(t) = 
F 2 (t ) = 0, 
E = 

G = 0 
H x = 
H 2 = 0 
H 3 = 


S 0 
0 5 


X 0 

0 x 


0 —fJL 

—p 0 


0 0 
0 0 


(43) 


We fix the uncertainty level at /x = 0.1. We then solve the associated //^ estimation problem 
using our Riccati equation approach. We again choose the desired disturbance attenuation level to 
be 7 = 0.65. The scaling parameters are chosen to be ei = 0.7 and e 2 = 1. A robust HA estimator 
is obtained as in (23) with the following parameters: 




0.3231 -1.3660 

-1.3660 0.3231 


B k 


— 1.4852t 
1.4852*. ’ 


C K = { 0.1 -0.1 ] . 


(44) 


The transfer function (33) of the estimator is obtained to be: 


—0.297ts — 0.3098* 
K ^ S ’ ~ s 2 - 0.6461s- 1.762" 


(45) 


The disturbance to error transfer function (35) of the robust filter in this example is shown in 
Fig. 3. We also plot those of the optimal H ca filter and the Kalman filter for the uncertain system 
(14) with d = l here for comparison. The optimal and Kalman filters are constructed for the 
nominal system of (14). 

Clearly, the robust Hoo filter again provides for preferable disturbance attenuation and 
performance compared with the standard optimal Hoc filter or the Kalman filter for this value of 
the uncertain parameter <5. 
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Bode Diagram 



7. CONCLUSION 

This paper considered the problem of robust estimation for uncertain linear quantum systems. 
The estimator is a classical filter, that produces a classical estimate of a variable of the quantum 
plant. The estimation problem is solved by converting it to a scaled IL X , control problem. The 
solution is obtained in the form of two algebraic complex Riccati equations. We have illustrated 
the results obtained by means of some numerical examples involving dynamic optical squeezers. 
As part of future work, the robust estimator constmcted here could be applied in studying 
robust coherent-classical estimation. For example, it might be interesting to explore if and when a 
robust Hoo coherent-classical estimator (with and/or without coherent feedback) can provide better 
estimation precision than the robust purely-classical estimator considered in this paper. Such a 
comparison of optimal estimators was presented in Ref. [32]. 


ACKNOWLEDGEMENT 

The first author would like to thank Dr. Obaid Ur Rehman for useful discussion related to this work. 


REFERENCES 

1. Wiseman HM, Milburn GI. Quantum Measurement and Control. Cambridge University Press, 2010. 

2. Yanagisawa M, Kimura H. Transfer function approach to quantum control-part I: Dynamics of quantum feedback 
systems. IEEE Transactions on Automatic Control 2003; 48( 12):2107—2120. 

3. Yanagisawa M, Kimura H. Transfer function approach to quantum control-part II: Control concepts and 
applications. IEEE Transactions on Automatic Control 2003; 48( 12):2121—2132. 

4. Yamamoto N. Robust observer for uncertain linear quantum systems. Physical Review A 2006; 74:032 107. 

5. lames MR, Nurdin HI, Petersen IR. H°° control of linear quantum stochastic systems. IEEE Transactions on 
Automatic Control 2008; 53(8): 1787-1803. 

6. Nurdin HI, James MR, Petersen IR. Coherent quantum LQG control. Automatica 2009; 45(8):1837-1846. 

7. Gough J, Gohm R, Yanagisawa M. Linear quantum feedback networks. Physical Review A 2008; 78:062 104. 

8. Maalouf AI, Petersen IR. Bounded real properties for a class of linear complex quantum systems. IEEE Transactions 
on Automatic Control 2011; 56(4):786-801. 

9. Maalouf AI, Petersen IR. Coherent H°° control for a class of linear complex quantum systems. IEEE Transactions 
on Automatic Control 2011; 56(2):309-319. 

10. Yamamoto N, Nurdin HI, James MR, Petersen IR. Avoiding entanglement sudden-death via feedback control in a 
quantum network. Physical Review A 2008; 78(4):042 339. 

11. Gough J, James MR. The series product and its application to quantum feedforward and feedback networks. IEEE 
Transactions on Automatic Control 2009; 54( 11):2530-2544. 


Copyright © 2016 John Wiley & Sons, Ltd. 
Prepared using rncauth.cls 


Int. J. Robust. Nonlinear Control (2016) 
DOI: 10.1002/mc 










14 


S. ROY AND I. R. PETERSEN 


12. Gough JE, James MR, Nurdin HI. Squeezing components in linear quantum feedback networks. Physical Review A 
2010; 81:023 804. 

13. Petersen IR. Quantum linear systems theory. Proceedings of the 19th International Symposium on Mathematical 
Theory of Networks and Systems, Budapest. Hungary, 2010; 2173-2184. 

14. Walls DF, Milburn GJ. Quantum Optics. Springer-Verlag: Berlin; New York, 1994. 

15. Gardiner C, Zoller P. Quantum Noise. Springer: Berlin, 2000. 

16. Wiseman HM. Doherty AC. Optimal unravellings for feedback control in linear quantum systems. Physical Review 
Letters 2005; 94:070405. 

17. Nurdin HI, James MR, Doherty AC. Network synthesis of linear dynamical quantum stochastic systems. SIAM 
Journal on Control and Optimization 2009; 48(4):2686—2718. 

18. Mabuchi H. Coherent-feedback quantum control with a dynamic compensator. Physical Review A 2008; 78:032 323. 

19. Sarma G, Silberfarb A, Mabuchi H. Quantum stochastic calculus approach to modeling double-pass atom-field 
coupling. Physical Review A 2008; 78:025 801. 

20. Petersen IR. Notes on coherent feedback control for linear quantum systems. Proceedings of the Australian Control 
Conference, Perth, Australia, 2013; 319-324. 

21. Bachor H, Ralph T. A Guide to Experiments in Quantum Optics. 2 edn., Wiley-VCH: Weinheim, Germany, 2004. 

22. Wiseman HM, Milburn GJ. All-optical versus electro-optical quantum-limited feedback. Physical Review A 1994; 
49(5):4110-4125. 

23. Lloyd S. Coherent quantum feedback. Physical Review A 2000; 62:022 108. 

24. Gough JE, Wildfeuer S. Enhancement of field squeezing using coherent feedback. Physical Review A 2009; 
80:042 107. 

25. Petersen IR. Coherent-classical estimation for quantum linear systems. Proceedings of the Australian Control 
Conference, Perth, Australia, 2013; 371-376. 

26. Roy S, Petersen IR, Huntington EH. Coherent-classical estimation versus purely-classical estimation for linear 
quantum systems. Proceedings of the Conference on Decision and Control, Los Angeles CA, USA, 2014; 5782- 
5787. 

27. Miao Z, James MR. Quantum observer for linear quantum stochastic systems. Proceedings of the 51st IEEE 
Conference on Decision and Control, Maui, Hawaii, 2012; 1680-1684. 

28. Fu M. DeSouza CE, Xie L. Hoo estimation for uncertain systems. International Journal of Robust and Nonlinear 
Control 1992; 2:87-105. 

29. Bouten L, Handel RV, James MR. An introduction to quantum filtering. SIAM Journal on Control and Optimization 
2007; 46(6): 2199-2241. 

30. Petersen IR, Anderson BDO, Jonckheere EA. A first principles solution to the non-singular H°° control problem. 
International Journal of Robust and Nonlinear Control 1991; 1:171-185. 

31. Zhou K, Doyle JC, Glover K. Robust and Optimal Control. Prentice-Hall: New Jersey, USA, 1996. 

32. Roy S, Petersen IR, Huntington EH. Coherent-classical estimation for linear quantum systems. ArXiv: 1502.03729. 


Copyright ©2016 John Wiley & Sons, Ltd. 
Prepared using rncauth.cls 


Int. J. Robust. Nonlinear Control (2016) 
DOI: 10.1002/mc 


